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The paper shows how a generalization of the elasticity theory to four dimensions and to space- 
time allows for a consistent description of the homogeneous and isotropic universe, including the 
accelerated expansion. The analogy is manifested by the inclusion in the traditional Lagrangian of 
general relativity of an additional term accounting for the strain induced in the manifold (i.e. in 
space-time) by the curvature, be it induced by the presence of a texture defect or by a matter/energy 
distribution. The additional term is sufficient to account for various observed features of the 
universe and to give a simple interpretation for the so called dark energy. Then, we show how the 
same approach can be adopted back in three dimensions to obtain the equilibrium configuration 
of a given solid subject to strain induced by defects or applied forces. Finally, it is shown how 
concepts coming from the familiar elasticity theory can inspire new approaches to cosmology and in 
return how methods appropriated to General Relativity can be applied back to classical problems 
of elastic deformations in three dimensions. 
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1 Introduction 

In the current perception of the general public (and of most physicists) nothing is as far apart 
as cosmology, on one side, and everyday life and its related physics, on the other. Of course the 
basic laws are the same, but the scales are so different that, when considering the universe at 
large, properties and phenomena emerge which are completely irrelevant at the human scale. It is 
however true that, when trying to describe and interpret what we see (or our instruments see) in 
the sky and deep in the past, we necessarily use the tools of reason with which we are equipped as 
well as the experience we draw from the world around us. Here we shall show how a phenomenon 
as familiar as elasticity and its related theory can be useful in the study of cosmology, and in return 
we shall apply a technique used in general relativity to the solution of simple ordinary problems 
of elastic deformation. The purpose is to stress the unity of physical sciences and the fruitfulness 
of the cross-fertilization among different areas of scientific and practical knowledge. Our start will 
be the observation that space-time, which is one of the basic ingredients of general relativity, looks 
like a four-dimensional physical continuum. By physical continuum we mean something endowed 
with peculiar properties and able to interact with the other ingredient which is matter/energy. The 
configuration of space-time and the way it reacts to the presence of matter/energy distributions are 
described in terms of geometrical properties. The same can be said of ordinary material continua 
in three dimensions. In both cases the best description of the properties and behaviour of our 
manifolds is cast in the form of tensors. The obvious difference is that for ordinary substances in 
three dimensions we have a dynamics expressed in terms of an independent parameter, time, which 
is instead part of the manifold when considering space-time. If we want to maintain a parallelism 
between material continua and space-time we shall compare 'equilibrium' configurations for both 
cases, i.e. 'static' configurations which may be labeled as being of equilibrium, respectively in three 
and in four dimensions. Now, how can we distinguish an equilibrium configuration from any other 
situation? A universally used technique is to write down the action integral appropriated to the 
system under study. Equivalcntly, we must consider the Lagrangian density of our manifolds. In 
fact, the approach we are describing is based on Hamilton's principle, which is a sort of 'economy' 
criterion not directly descending from other basic laws nor emerging from experiments. However, 
it is elegant and, above all, it works both for classical and quantum fields. In the case of General 
Relativity (GR), the Lagrangian density is chosen to be the simplest scalar which can be built 
from the geometrical tensor expressing the global and local configuration of space-time, i.e. from 
the Riemann tensor. That quantity is the scalar curvature R. In GR the action for an empty 
space-time is the Einstcin-Hilbert action, namely: 



four dimensions. As for R, its meaning is better understood when we remark that for an ordinary 
bi-dimensional surface it is proportional to the Gaussian curvature, i.e. the inverse product of the 
two principal curvature radii at a given point. The Lagrange coordinates of our system, i.e. the 
generalized coordinates that identify the state of the manifold (space-time), are the components of 
the metric tensor, i.e. the basic tensor representing the geometric properties of the manifold. The 
scalar curvature R is a function of the partial derivatives of the elements of the metric tensor up 
to (linearly) the second order. In a sense we may set a parallelism between the kinetic term in the 
classical three-dimensional Lagrangian (the kinetic energy), which contains the derivatives of the 
coordinates with respect to the external parameter (time), and R in four dimensions. However, if 
we imagine a physical continuum in three-dimensions obtained by deformation from an unstrained 
state we also have to take into account the elastic deformation energy. By analogy, extending this 
fact to four dimensions, we are led to adjust the Lagrangian density: to R we sum a scalar function 
of the Lagrangian coordinates accounting for the strain energy of the manifold. This approach is 
the essence of the Strained State Cosmology (SSC) theory (see Tartaglia [I], and Tartaglia and 
Radicella [5]). The explicit form of the additional energy density W will be worked out in the 
following. Once the new form of the action integral has been written, we may apply the stationary 
action principle and see what happens. We find a good description of an accelerated expanding 
universe: then, we shall verify that the same method works in three dimensions, also reproducing 
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known results usually obtained by other means. 



2 Lagrangian approach to cosmology including the strain of 
space- time 

In order to visualize the practical meaning of our approach, we may imagine two separate manifolds: 
one is a copy of the initial state of the actual manifold of our interest (reference manifold); the 
other is the natural manifold in its given deformed state. This picture implies that we adopt 
a higher dimensional point of view. For instance, if the manifolds we are studying are simple 
tridimensional surfaces, we look at them from within a three-dimensional space: sec FigfTJ where 
also the correspondence between pairs of positions on the two manifolds is shown. The embedding 
space is flat, which means that for this space Euclidean geometry holds; the same holds true for 
the reference manifold, but the natural manifold will in general be curved. With more dimensions 
(say N) the scheme is the same and the embedding manifold will be TV + n dimensional. If this 
description is used for space-time, the natural manifold, besides being four-dimensional, will include 
a Lorentzian signature, i.e. its tangent space at any location will be Minkowskian; putting it in 
simpler wording, the local speed of light will be the same everywhere, as required by relativity. 




Figure 1: Embedding space containing a flat and a curved manifolds 

We may now consider pairs of corresponding positions arbitrarily close to each other on the 
two manifolds. Let us compare the corresponding infinitesimal lengths. On the reference manifold 
we have 

d^ 2 = d^dC (2) 

The choice of the coordinates on the reference is arbitrary, since the result is a true scalar quantity, 
which means precisely that it does not depend on the coordinates. are the components of the 
Euclidean metric tensor in the representation appropriate for the chosen coordinates. 
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On the natural manifold the metric becomes non-Euclidean^]: 



d/' 2 = g„ v dx ti dx v (3) 

Now, are the components of the metric tensor on the natural manifold expressed in the coor- 
dinates chosen there. 

If we suppose that the natural manifold is obtained from the reference one by continuous 
deformation, we conclude that each position on one manifold has a unique correspondent on the 
other. If it is so, given an arbitrary choice of coordinates for one of the manifolds, an invcrtiblc 
function will give the corresponding values of the coordinates on the other. In practice, it is possible 
to express both line elements in terms of just one system of coordinates, for example the one used 
for the natural manifold. 

We can now express the difference 5 (dl 2 ) = d^ 2 — d/ 2 between the two line elements ^ and 
([3]) in the same coordinates: 

dr 2 - d^ 2 = (g^ - E^)dx"dx v = 2e^dx fl dx u (4) 
Formula defines the strain tensor e^ v . It is: 



(5) 



2.1 Einstein equations including the strain of space-time 

As we have already written, the intrinsic distortion of the natural manifold makes us wonder 
about the existence of some sort of deformation energy, which must be part of the Lagrangian as 
an additional potential (and in the equations it has to lead to a dynamical history of the universe as 
well). Starting from the ansatz that our space-time is homogeneous and isotropic]!, and assuming 
it behaves linearly, we may apply Hookc's law (whatever the number of dimensions is) in the form: 

= Miiv^ + 2/U£ MI , (6) 

where A and /i are the two Lame coefficients, accounting for the properties of the continuum, and 
a^v are the components of the stress tensor; rj^ are the components of the metric tensor of the 
flat space-time, i.e. the Minkowski space-time, which is the tangent manifold at any given event 
of the actual 'natural' manifold. 

The deformation energy density can then be written as 

W = (7) 
Casting © into (J7J the final form for the deformation energy is 

W= ^\e 2 + fie a pe af3 (8) 

In practice, we treat space-time as a physical continuum endowed with properties analogous to 
the ones of ordinary elastic solids. As a first step, we reconsider the Einstein-Hilbert Lagrangian 
density, C = R^/^^g, to which the deformation energy density must be added. Hence, the complete 
Lagrangian density in vacuo is 

C = {R + W)^—g (9) 



3 We shall remark that only the reference manifold is flat, thus the following relation holds only for the metric 
E^y (y indicates Cartesian coordinates): 

dy<* dy? 

" d «P 

and not for the metric g^v 

4 One remark about homogeneity and isotropy of space-time can be made. In fact, the unstrained manifold is 
obviously isotropic, since neither exogenous actions nor intrinsic defects are present. When we pass to the natural 
strained manifold, the locally anisotropic strain can induce some anisotropy in the elastic parameters. Nonetheless, 
this anisotropy may be considered as of second order with respect to the strain, thus the theory we consider is linear 
and the conditions of homogeneity and isotropy hold. 
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In general R 7^ 0, because we assume the presence, somewhere in the manifold, of a defect 
that generates intrinsic distortion. From (J5J) new generalized Einstein equations can be written 
differentiating with respect to the which are the Lagrangian coordinates of the configuration 
of space-time. The deformation energy W contributes to a new additional stress-energy tensor, 
T efi]y , which has to be read as an effective elastic stress-energy tensor. Hence, the final equations, 
the Einstein equations including the strain of space-time, look like 

Gnu = kT m „ + T efil , (10) 

Because T efiu is obtained varying a scalar with respect to a true tensor (the metric tensor g^ v ), it 
is also a good tensor, with all the properties of tensors. In explicit form: 

T e(iv = Xee^ + 2fj,e^ v (11) 

This equation comes from formula © after making the linear dependence of the stress from the 
strain explicit according to Hooke's law ([6]). Our elastic tensor is partially built from the metric, 
but it also constitutes an additional source together with the ordinary matter/energy term. In 
vacuo the Bianchi identities bring the Einstein tensor to zero implying the conservation of T eiMV , 
whereas in the presence of matter/energy the whole quantity kT^ + T efiV is conserved, signifying 
the possibility of a transfer of energy between the matter and the strain, which is a well-known 
subject in classical physics. 

2.2 A Priedmann-Lemaitre-Robertson- Walker universe 

The strategy described so far has been applied to the study of the properties of the universe as a 
whole, assuming a global Robertson- Walker symmetry, i.e. homogeneity and isotropy in space (over 
the scale of hundreds of megaparsecs (Mpc)). The result is visually summarized in Fig 121 where 
the bell-shaped surface represents the space-time of our Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe, including a defect ('initial' singularity) and displaying (horizontal cross-sections) 
an expanding space, whose expansion rate is initially decreasing, then accelerating. 

The flat surface on the top of the figure represents the Euclidean manifold from which, according 
to SSC, the real space-time has been obtained by deformation induced by the presence of a defect. 
The cosmic time is measured from the singularity (the defect) along the generators of the bell- 
shaped manifold (dashed line) . The circles and the arrows shown in the figure manifest the peculiar 
gauge choices corresponding to different deformation strategies. 

3 The curvature-plus-strain method applied to 
three-dimensional solids 

We are now demonstrating that the method outlined in the previous pages might be successfully 
applied to three-dimensional material continua. Obviously, the time dimension is off the table: for- 
getting about space-time, we remain with a space in three dimensions with no Lorcntz signature. 
However, the idea which we are trying to convey is that solids can be seen as three-dimensional 
manifolds that turn from their unstrained reference state to a curved natural configuration. We 
represent the unstrained state by means of a flat manifold, whose metric is determined hereafter. 
This is indeed a non-conventional interpretation of the theory of elasticity. The non-conventional 
aspect is not in the tensorial machinery, which was originally developed having continuum me- 
chanics in mind, but in the use of the Einstein-Hilbert Lagrangian density with the addition of 
the deformation energy density. We may look at our approach as at an exercise, in which we shed 
light on the relations existing between two disparate branches of physics: cosmology and continua 
mechanics. 

Hence, we are about to consider solid bodies with particularly symmetric shapes, as well as 
loaded symmetrically Notably, we will develop the case of central symmetry, i.e. a spherical cavity 
immersed in an infinite medium, where the stresses originate from uniform and isotropic loadings 
which obey the symmetry of the body. For this case the classical theory of elasticity provides exact 
solutions that we are using to prove our thesis. 
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Figure 2: A homogeneous and isotropic expanding universe, obtained by continuous deformation 
from a flat Euclidean manifold with a defect (the black disc on the reference manifold which maps to 
a single point on the natural manifold). The arrows correspond to different deformation strategies 



3.1 Metric tensors for the natural and reference manifolds 

The natural manifold is a three-dimensional material continuum, which in general is not flat, but 
it has assumed a curvature as a consequence of some action applied to it. The reference manifold is 
an ordinary three-dimensional flat space: with no loss of generality we assume it coincides with R 3 , 
described with a convenient set of coordinates, spherical in our case. Consequently, the embedding 
space has to be at least in four dimensions to contain the manifolds described above: the simplest 
possible is perhaps the four-dimensional Euclidean space R 4 . 

In the fashion of the SSC, when a solid is unstrained its metric is flat, and in spherical polar 
coordinates this means: 



Because the loading respects the central symmetry of the sphere, after deformation the symmetry is 
preserved, even though dimensions are slightly scaled. Fortunately, only one scale factor is needed 
to describe the metric, and it is function of the only coordinate r: we call it / = f(r). The freedom 
in choosing the coordinate set offers many possibilities for the explicit form of the line element, 
always preserving the symmetry. Here we choose: 



Consider now the reference manifold. It is flat, and its metric E^ v is the tensor needed to 
represent the length of a line element before deformation. The E )n , is indeed Euclidean: however, 
its explicit expression in terms of the coordinates used on the natural manifold involves knowing 
how the material was deformed. 

Figf2] gives a clear idea of the problem. Point P on the natural manifold may correspond to 
infinitely many points on the reference manifold. For example, assuming that P lay on the solid 
circle before any stress was applied, the material shrank more than if it were on the dashed circle, 
but less than if it were on the dash-dot one. The various ways to go from the strained solid 




(12) 




(13) 
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backward to the reference flat state depend on the strategy of deformation adopted. This strategy 
is represented by the function w = w(r), which in practice establishes a one-to-one correspondence 
between points on the two manifolds. The choice for the correspondence is in principle arbitrary, 
but it can be determined by looking for the less energy consuming strategy. The coloured arrows 
in Fig 12] are the effective vectors of a displacement field u: we may think of them as strings that 
pull points of the flat manifold, moving them to the natural configuration. Therefore, the same 
line element on the natural manifold may have several corresponding line elements on the reference 
manifold. Of course we must fix the form that w assumes, eventually. In conclusion, the function 
w simply defines the scale of the line element on the flat reference space, that has to be compared 
with the corresponding line element of the natural frame. Summing up, the metric E^ v of the 
reference manifold is 

ds 2 = w l2 {r)Ar 2 + w 2 {r) (dtf 2 + sin 2 § dip 2 ) (14) 
where ' stands for the total derivative with respect to r. 



3.2 Lagrangian density 

The Lagrangian density C is indeed similar to but with the more general \f\g~\ instead of ^J—g: 



C = (R + W) \f\g\- The two ingredients R and W are ready to be built, because the metrics are 
already well defined. 

The curvature of the reference manifold is obviously zero, whereas the curvature of the natural 
manifold is directly computed using the metric tensor g^ of (|T5j) : 

2rf + 2/ 2 -2f , . 

A= r 2f2 ~ (15) 

On the other hand, applying (JS|), the non-zero components of the strain tensor in covariant 
form are 

f - w' 2 

e rr = L__ (16) 

2 2 

e<H> = ^ (17) 



2 2 

£ipip = ^ sm2? 9 (18) 

To raise one subscript of the strain tensor in order to obtain the mixed form e" let us use the 
inverse metric g^ u : 

f - w' 2 

4 = (19) 

2 2 

r — w 



In this way, the scalar strain is immediately built: 

£ = £; = iz^! + r!^! (21) 

Applying again the inverse metric g^ v to e", we easily get . So, we can compute the second 
degree scalar 



(f-w' 2 ) (r 2 ~w 2 ) 



The deformation energy is then obtained by casting these results into ((SJ). Eventually, the deter- 
minant of the metric yields 

\ 7'-' -mi' (23) 
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And the Lagrangian density is: 



C 

sintf 




ylr 2 (24) 



4 Results: extremization of the Lagrangian 



From (|24[) . we could develop the so-called Lagrange equations. Nonetheless, these equations would 
be highly non-linear and very difficult to solve. There exists a simpler way to verify that a 
Lagrangian density like © accounts also for the deformation of a three-dimensional solid. Recall 
that, in the case of a field theory, Hamilton's principle can be stated as in Landau and Lifshitz [3] 
(pp.287-288): 



SS 



Cd 3 x 



fit 3 ) 







(25) 



Hamilton's principle (|25[) states that the action must be extremal (a minimum in general, but not 
always so): then, it also involves that in order for SS to be zero, the Lagrangian density C is 
necessarily extremal. This is why Hamilton's principle will be restated as 



SC = 



(26) 



which means that when we perform the variation of the Lagrangian density we find zero. 

Where w and / are being derived from? We can use the exact solutions worked out within 
the classical elasticity theory. In particular, we wish to determine the strain tensor for a hollow 
spherical cavity of radius b subjected to a spherically symmetric loading, say a uniform internal 
pressure p. The walls of the cavity are assumed to be infinitely thick, so that no stress acts from 
outside. We use the spherical polar coordinates (r,§,(p), and note that the displacement vector 
u is function of the only variable r. This implies V x u = 0, and the Navier-Cauchy equilibrium 
equation yields: V • u = constant. We integrate this equation, and use the strain-displacement and 
stress-strain relations (following the approach outlined in Landau and Lifshitz [4], pages 20-21). 
Imposing boundary conditions as follows 



we find that the strain tensor is 



e i9 



if r 
if r 



pb 3 1 
2/7 ^ 
pb 3 
"4^ 



1 



(27) 

(28) 
(29) 



Everything here is expressed using coordinates defined with respect to the unstrained reference 
frame. This is the common practice when studying elasticity; in the cosmological application 
instead, the natural choice is to use the coordinates defined for the actual universe where we live 
in, which correspond to the strained rather than the reference frame. In any case, as we saw before, 
for the final purpose the two choices are perfectly equivalent. Comparing (|28[) and (|29[) in covariant 
form with (ITBl and (fT7|) we get 



/ = 



w = 



2r 3 p + pb 3 



8pb 3 r 3 p, + 9p 2 b 6 



2pr 

- 16r 6 /z 2 



-pb 3 ) 



(30) 
(31) 
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In practice, we have rewritten the exact solutions in terms of / and w, which we now cast into 
thus finding an explicit expression for the Lagrangian density. Let us call it C. Performing 
the variation on C(r), we should find that 5C = holds up to higher order terms. This fact is 
apparent if we consider that the SSC is not linear, but the values for / and w we cast inside are 
obtained from a linear theory. 

In other words, what we are proving now is that these solutions minimize the Lagrangian 
density. To do so, we vary the values of / and w by perturbing the exact solutions ([25)1 and (|29l) . 
and find that the Lagrangian density immediately increases. For instance, consider the perturbation 
e rr +a(r), where a(r) is a small function (compared to s rr ), and is left unperturbed. Proceeding 
as before we have 

r 2 — w 2 pb 3 1 f — w' 2 pb 3 1 

^ = — ; an d = — o + otir) 

2 Afir 2 2fi r 3 w 

Solving / and w, and plugging these values into C(r) (expanding with respect to powers of a) we 
find the perturbed Lagrangian C: 

£ = £ + a -(0) + a 2 •(...) + ... 

Although up to first order C = C, if we add a further term in the expansion, i.e. second order terms 
(a 2 collects these terms, whereas a collects none), the Lagrangian does get bigger, indeed. These 
results are shown in FigGH where the computation has been repeated moving the perturbation onto 
the other non-zero strain component, i.e. + a(r). 

Our simple example with three-dimensional spherical symmetry corresponds in cosmology to 
a closed ever expanding FLRW universe. If we had treated a different symmetry, also the four- 
dimensional correspondence would have been to other cosmological models with different topologies 
of the initial singularity. 



5 Conclusions 

Wc have shown how a generalization of General Relativity, obtained exploiting the analogy with the 
ordinary linear theory of elasticity in three dimensions, is able to account for the global behaviour 
of the universe, and in particular to reproduce the accelerated expansion discovered in 1998 and 
1999 observing the luminosity of type la supernovae (Riess et al. [5], and Pcrlmuttcr ct al. [5]). 
The situation is visually described in Figj2] where successive spacelike cross-sections of space-time 
increase at a growing rate with the distance from the initial singularity (i.e. with cosmic time). 
This result is contained in ref.s [3] and [7|. 

Our result is obtained adding to the classical Einstein-Hilbert action for empty space-time a 
term accounting for the equivalent of a deformation energy, built on the strain induced in space- 
time by either the presence of texture defects (as some times happens in ordinary three-dimensional 
continua) or of matter/energy distributions. The additional term lends itself also to a simple 
interpretation of the so called 'dark energy' of cosmology. Working out formally the equation of 
state of the strain term as if it were a fluid in the FLRW universe, we see that the w parameter is 
negative (for big enough values of the scale parameter a) , which corresponds to a negative pressure 
(ref.[2], sec. 5): a property which is ascribed to the dark energy. More precisely we should write 
that the strain term resembles the 'quintessence' of cosmology, since the parameter of state is not 
equal to -1 (cosmological constant) but corresponds to a function of the scale factor. The fact that 
we refer to a strain fluid rather than to a solid is not important here, because the solutions we find 
in four dimensions are the equivalent of equilibrium configurations so that the distinction is only in 
the explicit form of the energy density. The difference between the classical picture of an evolution 
in time and the four-dimensional manifold of GR is in that an evolution scenario would require 
at least a fifth dimension with an evolution parameter replacing the role of the independent time 
of classical dynamics. This is why we say that any global solution for the configuration of space- 
time is an 'equilibrium' because it does not depend on any further evolution parameter, just as it 
happens for the equilibria of three-dimensional classical systems. The strained state is represented 
by the strain tensor, intended, as it is the case in the theory of elasticity, as half the difference 
between the intrinsic metric tensor of the actual deformed continuum and the metric tensor of the 
initial unaffected manifold. Whenever the strain is identically zero we have a flat manifold, i.e. a 
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Exact Lagrangian Density and its Perturbations 



Exact and order 1 
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Figure 3: Lagrangian density as a function of the radius (dash-dot line) and perturbations thereof 
expanded up to second order (the dashed curve represents the case where e rr is perturbed, and 
the solid one the case where is perturbed). When we expand the perturbed Lagrangians only 
up to first order, the curves coincide with the dash-dot one 

continuum whose intrinsic geometry is Euclidean. In the case of space-time, gravitational effects 
are included in the strain tensor, i.e. in the non-trivial part of the metric tensor. Once the action 
integral has been generalized in the way outlined above, the rest is deduced applying Hamilton's 
principle. Next, we have shown that the same approach adopted at the cosmological scale is viable 
also in three dimensions to treat ordinary problems of elasticity. Here too, even though more direct 
approaches are often at hands, the equilibrium configurations of elastic continua may be deduced 
from the three-dimensional version of the SSC action integral used in cosmology. The reference 
to cosmology is rather appealing, but we would like to stress that the method we have outlined 
is indeed an extension of GR and can be applied to typical problems of gravitation with different 
symmetries. It has for instance been applied recently to cylindrical symmetry in four dimensions, 
i.e. to the Schwarzschild problem (spherical symmetry in space) [5] and will be applied to screw 
symmetry (rotating systems). We hope that the corresponding lower dimensional symmetries can 
help in solving the GR problems. 

Summing up, we have seen an example of how the universal nature of physical laws is manifested 
establishing a connection between the universe as a whole, or at least space-time that is the arena 
on which the cosmic drama is played, and the everyday familiar description of the properties of 
deformable solids. 
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